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Abstract: Recent calculations in one-loop and Gaussian approximation, using the so-called 
autonomous renormalization scheme, indicate a comparatively massive, narrow Higgs excitation 
at about 2 TeV. Here I show that this result qualitatively persists in the framework of a post- 
Gaussian variational approximation for the pure 0(iV)-symmetric 4 -theory for N > 1. The 
method is based on nonlinear transformations of path- integral variables, and the optimization 
amounts to a Schwinger-Dyson-type summation of diagrams. In the case of 0(4), for example, I 
find Afjjiggs = 2.3 TeV, compared to 1.9 TeV and 2.1 TeV in 1-loop and Gaussian approximation, 
respectively. 
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1 Introduction 

Five years ago Stevenson and Tarrach [jlj] suggested certain UV-nows for the bare parameters of 
the </> 4 -model which lead to a finite Gaussian effective potential (GEP) H in 3 + 1-dimensions. 
The method was termed autonomous renormalization scheme (ARS) because it is not directly 
related to the perturbative renormalization, the essential novelty being an infinite rescaling of 
the field. The autonomously renormalized GEP exhibits all generic features necessary for the 
Higgs mechanism, but it eludes triviality in that it is not equivalent to the tree-level potential. 

Recently essential progress has been made in understanding the theoretical basis of the 
ARS and in deriving predictions for physical quantities from it. In a series of papers Consoli 
et al. j|] investigated the ARS with renormalization group methods. They found asymptotic 
freedom in the asymmetric phase. Further, they discovered that the ARS also works in the 
loop approximation on the one-loop level and that definit predictions for the Higgs mass can be 
obtained under the natural assumption that the theory has a massless particle in the symmetric 
phase ||. 

Compared with standard methods, the ARS predicts a relatively high Higgs mass at about 

2 TeV, varying in the 10%-range, depending on whether the 1-loop potential or the GEP is 
considered, or whether couplings of the scalars to gauge bosons are taken into account or not 
|H|. For instance for the 0(4)-theory in Gaussian approximation one finds MHi gg s= 2.05 TeV. 

One of the most important questions which must be asked in this context is whether the 
autonomous ARS has a realization beyond GEP and one-loop approximation. This is precisely 
the point which is tackled in this paper with a post-Gaussian variational calculation. During 
recent years a method has been developed, which allows to generalize the variational approach 
in quantum field theory to non-Gaussian trial states in the canonical formalism (?], || and non- 
Gaussian trial actions in the path-integral formalism Q, respectively. In the following I shall 
work in the covariant setting and restrict myself to the pure 0(7V)-theory. 

2 The method 

The generic (euclidean) action considered throughout the paper is given by 

S = ~ J (-dvpdvtf + m|0V) + A B J ( W 2 (1) 

with J = J d 4 x. In order to approximate the effective potential, V[0 C ], I introduce the nonlinear 
transformation 

1 = fj l (p) + Xo S( p ) + s c(q,r)fi a (q)f] a (r)5(p-q-r) 

Jqr 

4> a = fj a (p) a = 2,...,N (2) 



which has Jacobian equal to unity, where / = l/(2ir) 4 J d 4 p, S(p) = (2tt) 4 5(p), and tildes 
indicate Fourier amplitudes. 

Transformation (|2|) is appropriate for studying spontaneous symmetry breaking O(N) — > 
0(N — 1) and has previously been used in the canonical approach |g|. Leaving the path-integral 
measure invariant, (g) is well suited for a variational calculation, because it can be evaluated 
in closed form, i.e. does not lead to a series expansion of the expectation value. The c-number 
Xo an d the correlation function s c(q, r) - the factor s is split off for normalization and as a 
bookkeeping device - are variational parameters. 

The upper bound on the effective potential derived from (ffl) is 

V \<j> c ) <V A {4>c) = -^r mm j- log Af+AT 1 j D n e~ SG ^ (S[<f\ - S G [r,})^ , (3) 

where the available variational parameters have to be optimized under the constraint 

(f> c := AT 1 f Dr] e~ Salri] <f) 1 = fixed, (4) 

i.e., with the expectation value of the field held fixed, and 

AT = f D V e- Sa[r > ] . (5) 

In the equations above the </>'s have to be substituted with transformation (S). Sg is a quadratic 
test action, given by 

SgM = \ J {vHp) Gl\p) fj\-p) + fj a (p) G T Hp) f, a (-p)} , (6) 

where Gl and Gt are adjustable propagators for (with respect to the direction of symmetry 
breaking) longitudinal and transversal modes, respectively. 

3 Autonomous Gaussian approximation 

If the nonlinear term is absent in (0), the simple identy 4> c = xo holds and the optimization 
yields the (covariant) GEP. (In contrast to the canonical formalism, this procedure can be easily 
extended to the effective action ||, but for my present purpose the effective potential is suffi- 
cient.) In the following I am summarizing the basics of the ARS for the Gaussian approximation, 
and afterwards the analysis is extended to the nonlinear transformations. 
The bare unoptimized GEP is given by 

V G = J L + (N - 1) J T + im| [l L + (N - 1)I T + (g 

+A B [3(I l + <fcf - 2cpj + 2(N - 1)I L (I T + 4>l) + (N 2 - 1)(/ T ) 2 ] , (7) 

with 



J L ^ = \J p lo g G- L \ T) + \J p P 2 G L{T) (p) 



(8) 



and 
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(10) 

p~ -|- HIT 

with optimization equations 

Q 2 = ml + 4A B (3J (n) + 3^ + (N - l)I (u>)) (11) 

for longitudinal and 



ui = m 



| + 4A B (/o(^) + ^ + (A r + l)/oM) (12) 



transversal mass M. With (10) the J's and /'s become 



J(m) = -- log G (p;m) - -m 2 I (m) , I (m) = G (p;m) . (13) 

The ARS consists of the following UV- flows for the bare parameters: 

^-12aIo(0) ) ^ = a ^ ^ = a 

J-i(/*) ^-i(w 

where 

/_ 1 ( / i) = 2/(G (p;^)) 2 (15) 

is a logarithmically divergent quantitiy which plays an inportant role in the following, and fi is 
an arbitrary renormalization mass, a dimensional transmutation parameter for it replaces the 
dimensionless coupling constant. The second parameter, niQ, which turns out to be the mass at 
the at the origin, is set to zero in the following ||. 
Using (|lj) in (j^) one obtains 

V G = al-i(fi) (Aft 4 + Blu 4 + CfiV + DQ 2 z <S> 2 + Eu 2 z <5> 2 + zg$ 4 ) + finite + £>(1/I_i) (16) 

with 

A = l + l B = ^ + ^T i - ° = W-V, D = -3, B=-(JV-1), (17) 

i.e. all quartic and quadratic divergencies are cancelled, leaving behind a logarithmically di- 
vergent expression. In order to end up with a finite effective potential eventually, also the 
logarithmically divergent term, proportional to I_i, has to vanish in the minimum with respect 
to the variational masses. This amounts to an equation for a which can be solved analytically. 
The result is given by: 

a = — \ . (18) 

4(1 + VFT3) 



The remaining finite terms in (0) are 

Vfc(*) = ^z 2 ^ (*log(<W) - Y) (19) 

with 

X = ^^ - (20) 

{2 + VWTs)(i + VWT3) 2 

For X > the potential is stable, has massless excitations at $ = and an asymmetric vacuum 
with massive particles. (The constant Y can be calculated but is not important for the following 
considerations.) The longitudinal mass - which is to be with the Higgs mass in the asymmetric 
minimum - is given by 

n 2 = z a L $ 2 with a L = 8a . (21) 

It is worthwile to mention that the stability condition X > is automatically satisfied, once 
the solution for a exists. In ( |l6|) all quantum corrections have to add up in order to produce 
a negative contribution that cancels the positive ZqQ 4 . The term proportional to X is solely 
generated by those subleading terms, which are produced when the longitudinal and transversal 
masses are eliminated in favor of /a. The latter is done with the help of the relation 

I Tfl 

/_i (m) = 7 n (/j) - — log — , (22) 



and the log-term yields the log(<3? 2 //j 2 ) in (|l9|). All other subleading terms contribute only to 
Y. Thus, X > is a direct consequence of the minus sign in (E2|). 



The finite factor zq that enters the wavefunction renormalization in (14) can be calculated 
from the renormalization condition 

d 2 V 



= n 2 



*=*„ = «*w ■ ( 23 ) 



d$2 

where & v is the value of the field in the asymmetric minimum. This leads to 



2 °L 2 10 + 6VF+3 + 2iV 

Zq = 7T — - = 7T . . (24) 



If now the well-known value <&„ = 0.246 TeV is inserted in (|2lD, Mniggs can De computed; the 
table below shows the results. 

4 Non-Gaussian contributions 

I am switching on the nonlinear term now, leaving the correlation function c(p, q) subject to 
optimization. The expectation value can be calculated straightforwardly. The unoptimized 
effective potential becomes 

V A = V G + AV = V G + s 2 X 7+2 m BS 2 X2 + X B {4<J) c sx3 + 2s 2 X5 (25) 

+s 2 [67 L + 2(N - 1)7 T + 6^] X2 + 4cP c s 3 X4 + 3s\ X6 + xl)} (26) 



where Vq is the Gaussian result, with yet undetermined longitudinal and transversal propagators, 
and 



X2 

Ipq 

X3 

Ipq 



2{N-l) f c 2 (p,q)G T (p)G T (q) 1 (27) 

Jpq 

2(N-1) f c(p,q)G T (p)G T (q), (28) 

Jpq 

X a = 8(N-1)[ c(p,q)c(p,r)c(q,-r)G T (p)G T (q)G T (r), (29) 

Jpqr 

X5 = 8(N-1)[ c(p,q)c(q,r)G T (p)G T (q)G T (r), (30) 

Jpqr 

X6 = 16(JV-1)/ c(p,q)c(p,r)c(q, S )c(r,s)G T (p)G T (q)G T (r)G T (s), (31) 

Jpqrs 

XI = (N-l) I {p + qfc 2 {p,q)G T {p)G T {q). (32) 

Jpq 

It is the term linear in s, proportional to X3i that assures an improved energy compared to 
the GEP, at least on the level of the bare theory. Neglecting for the present the higher-order 
terms proportional to X4 and Xe - the justification for that will be given below -, the optimization 
can be carried out, reducing the problem to a Schwinger-Dyson-type integral equation for the 
optimal correlation function c(p, q) : 

c(p, q) = G (p + q- Q) (l - 8A B J [c(p, r) + c(q, r)\ G T (r)^j , (33) 

with 

U 2 = m 2 B + 4A B ((N - 1)I T + 3(/ L + s 2 X2 + eg)) (34) 

and where 

s = -4\ B <Pc (35) 

has been chosen for convenient normalization. 

The solution of ( |33[) leads to a net energy-decreasing term 



AV = -8X1^X3 + 0(s 4 ) (36) 

where X3 is ( HI) with optimal correlator. The ^-integrals satisfy 

2x7 + ^ 2 X2 + 4A B X5 = X3 • (37) 

The iteration of ( |33| ) reveals that (36) contains an infinite series of contributions to the two- 
point vertex function - the first four graphs are depicted in Fig. 1 -, thus demonstrating that 
the approximation goes far beyond the Gaussian approximation which sums up the "cactus" 
graphs. 



5 Renormalization 

Dimensional analysis shows that the bare correction AV is quadratically divergent. As a conse- 
quence, it necessitates a correction to the Gaussian mass counterterm in (pi). As discussed in 



[13|, however, any such correction to m 2 ^ multiplies the whole mass term, 

I L + s 2 X2 + <g + (N - 1)I T ) , (38) 



thus generating new (quartically and quadratically) divergent terms which have no direct coun- 
terpart in the expectation value. This dilemma originally led to the notion of the instability 
of variational calculations [jllj] in general and the ARS in particular [|l^] . With the help of the 
diagrammatic representation, however, a deeper understanding of the problem and, eventually, 
a straightforward solution can be obtained. 

The first diagram in the series Fig.l, the so-called barred circle, has the close relative shown 
in Fig. 2, a vacuum diagram that may be generated from the barred circle by attaching two ex- 



ternal lines and connecting them. In second-order ^-expansion [1C], for example, both diagrams 
are present and the renormalization works as Fig. 2 contains the barred circle as a divergent 
subgraph. In the variational calculation, however, the situation is quite different. Fig. 2 is 
not generated automatically. This can be verified by considering the effective potential at the 
origin where the optimization can be carried out exactly (8|. Consequently, the counterterm for 
the renormalization of the barred circle causes (quartically) divergent subtractions which are 
actually supposed to renormalize vacuum graphs. 

There are two possible solutions to this problem. Firstly, one can think of an ansatz which 
generates both diagrams. This can in principle be achieved by a transformation like (J2|) with 
cubic nonlinear term of the form r\ rfrf. However, this transformation has nontrivial Jacobian 
and thus makes the evaluation of the ansatz much more complicated and the applicability of the 
variational method questionable. Secondly, as suggested in O], one can choose a sub-optimal 
form of the variational parameters, just generating the divergence, which otherwise would be 
caused by Fig. 2 and which then can be subtracted by the mass counterterm. 

The ansatz is given by 

G L {p) = G Q (p-n)-s 2 i{p) (39) 

with 

i{p) = 2(N-l)G T (p) f c\p,q)G T (q) (40) 

Jq 

for the longitudinal propagator and 

Gt(p) = G (p;oo) (41) 

for the transversal propagator. 



With (39) one obtains 



I L = Jo(O) - s 2 X 2 (42) 



and 



J L = J(Q) + Ul 2 s\ 2 + 0(s A ) (43) 



The next step is to set 

s = -4A B Uc ~ yf-A(n)-(N-l)A(u;)\ (44) 

where 

A(m) = I (m) - J o (0) = -\m 2 I^) + ^ hog ^- - 1 J . (45) 

As opposed to (^) the ansatz (^J) generates a term proportional to 

(A(n) + (JV-l)A(u;))x3, (46) 



whereas on account of the identity ( J37| ) terms containing the square root of the A's cancel among 
each other. 

Finally, by demanding that the excitation at the origin remains massless, the bare mass 

ml = -12A B /o(0) + 16A| X3^ c=0 (47) 

is obtained. As a result of fl39|), ([44]), and ( |47| ) the effective potential is given by 

v* = vfc - sa| (xs - xalo) (A(n) + (jv - i)A( w ) + ^) + 0(0*) . (48) 

This expression is not yet finite, but all power divergencies have been removed. As I shall 
demonstrate below, the only remaining divergencies are logarithmic and, like in the GEP-case, 
they can be cancelled by adjusting the parameter a in Ab- 

Before I proceed, the terms O(cj)f) have to be discussed. From ( |i~4|) and (|44]) it follows that 
s 2 ex l/i_i. Analyzing the integrals which occur in 0(cp A ) (like \ii X6> an d other contributions 
produced by (|39D), one finds that all these terms are at most (9(1/I_i) and, thus, they are not 
affecting the finite parts of the effective potential in the ARS. 

6 Evaluation of the effective potential 

In this section the Schwinger-Dyson equation (^) will be discussed. The equation is in a sense 
the higher-order analogue to the self-consistency equation of the Gaussian approximation. The 
integral equation corresponding to the latter is extremely simple and can be solved algebraically. 
With (|33| ) the situation is more complicated. The technical problems are closely related to the 
ones occuring in the \/N expansion |14||. Like there, subsets of of diagrams, containing chains of 
bubbles of different length, can in principle be summed up. In the following I am not proceeding 
along those lines, however, because there are good reasons to believe that the series converges 



rapidely when the factor a (which is the expansion parameter in the integral equation) assumes 
its physical value, where logarithmic divergencies are cancelled. 

The leading (logarithmic) divergence of the diagrams generated by (33) is proportional to 



I_i , where k is the number of vertices. With the help of (45) one finds that AV oc J_i, like in 
the Gaussian approximation. I have calculated explicitely the contributions from the first and 
second diagram in Fig. 1. The leading divergence of these graphs is given by 

- - (ft 2 + 2u?)l 2 _ 1 and - ~ (ft 2 + 2u; 2 )I-l > (49) 

8 3 

respectively. 

With the ansatz 

ft 2 = o L z § 2 , J 2 = a T z $ 2 (50) 

the divergent part of the effective potential again has the form ([!(]) . The coefficients read 

A' = A -A, B' = B - 2{N - 1)A, 

C' = C-(N + 1)A, D' = D + 4A E' = E + 2A (51) 



with 

A = (N - l)a (l - — a + 0(a 2 )) (52) 

and where the unprimed letters stand for the Gaussian coefficients (|l7]). One derives three 
equations from the conditions that the term proportional J_i vanishes in the minimum with 
respect to the variational masses. These equations can be solved numerically giving a, ctl, and 
a? as functions of N. 

Eventually also the finite parts of the effective potential can be obtained. They again have 
the form (|l9|), and the factor X can be calculated on account of (^) directly from the leading 
divergence. The results for the Higgs mass are given in the table. 

7 Discussion of results 

I have demonstrated that within the 0(N)- theory the autonomous renormalization scheme 
exists beyond the one-loop and Gaussian variational approximation. This result gives some 
confidence in the ARS and makes it unlikely that one is dealing with an artifact of the GEP. The 
optimized nonlinear transformations lead to a Schwinger-Dyson-type integral equation, formally 
summing up a series of diagrams that goes far beyond the one that is taken into account by 
the Gaussian method. The integral equation has been approximated by the second iteration, 
which contains diagrams with up to three loops. An estimate of the changes due to four-loops 
indicates that the error of the approximation is less than 3%. 

The results for the Higgs mass are shown in the table. For low values of N, the mass is 
higher than in one-loop and Gaussian approximation. For instance, for A^ = 4 it lies 10% above 

9 



the Gaussian and 20% above the one-loop prediction. For higher N it drops substantially below 
the Gaussian results. 



N 


Gaussian 


Non-Gaussian 


1 


2.19 


- 


2 


2.13 


2.21 


4 


2.05 


2.27 


7 


1.97 


2.22 


10 


1.95 


2.08 


20 


1.84 


1.70 


100 


1.69 


1.34 



Table 1 Values of Higgs mass (in TeV) for different values of N from Gaussian and non-Gaussian 
approximation . 



The present calculation shows that the error in the Higgs mass due to the approximative 
treatment of the self interaction lies in the same order of magnitude as corrections due to 
interactions with gauge bosons and fermions. Hence, if the phenomenological consequences of 
the ARS are to be taken seriously, one clearly has to confirm them by alternative improved 
approximations, like the conventional loop expansion, the (^-expansion fK|, or the optimized 



expansion 15]. 
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Figure captions 



Fig. 1 Four contributions in the series of diagrams generated by the Schwinger-Dyson equation 
). Thin (thick) lines represent transversal (longitudinal) propagators. 



Fig. 2 Vacuum diagram that contains the "barred circle" as divergent subgraph. 
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